Thermodynamic formalism for field driven Lorentz gases 
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We analytically determine the dynamical properties of two dimensional field driven Lorentz gases within the 
thermodynamic formalism. For dilute gases subjected to an iso-kinetic thermostat, we calculate the topological 
pressure as a function of a temperature-like parameter f3 up to second order in the strength of the applied field. 
The Kolmogorov-Sinai entropy and the topological entropy can be extracted from a dynamical entropy defined 
as a Legendre transform of the topological pressure. Our calculations of the Kolmogorov-Sinai entropy exactly 
agree with previous calculations based on a Lorentz-Boltzmann equation approach. We give analytic results for 
the topological entropy and calculate the dimension spectrum from the dynamical entropy function. 
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I. INTRODUCTION 

In dynamical system theory the Lorentz gas acts as a 
paradigm allowing to address fundamental issues of non- 
equilibrium processes. Recently, dynamical quantities of this 
system such as Lyapunov exponents or Kolmogorov-Sinai en- 
tropies in a non-equilibrium steady state have been calculated 
analytically |1, 2]. Systems like the Lorentz gas in a non- 
equilibrium steady state can be modeled on the microscopic 
level by introducing a so-called thermostat which removes the 
dissipated heat from the system. Assuming that fluids can be 
regarded as hyperbolic systems on a microscopic level one can 
relate the dynamical quantities to viscosities or diffusion coef- 
ficients OIIj,!^]- This connection is based upon a relation be- 
tween phase space contraction rates and entropy production. 

In the 1960s and 70s Sinai, Ruelle, and Bowen have de- 
veloped a formalism for dynamical system theory which due 
to its striking similarity to Gibbs ensemble theory was given 
the name thermodynamic formalism |6, 7,_8J. This formal- 
ism is applied to hyperbolic systems and like in ordinary sta- 
tistical physics a partition function is defined which is con- 
structed by giving points in phase space a particular weight. 
From this a central quantity, the topological pressure, is de- 
rived which is the dynamical system analog of the Helmholtz 
free energy. From the dynamical partition function properties 
like the Kolmogorov-Sinai (KS) entropy or the topological en- 
tropy are obtained through the pressure and its derivative with 
respect to a temperature-like parameter (3. Also derivable are 
dimension and entropy spectra and, for systems with escape, 
escape rates. 

The paper is organized as follows. In Sec. |n]we briefly re- 
capitulate some properties of the thermodynamic formalism. 
Sec. Hn] introduces the field driven Lorentz gas and Sec. IIVI 
the concept of the radius of curvature. Calculations for the 
field driven random Lorentz gas within the thermodynamic 
formalism in two spacial dimensions are presented in Sec.lVl 
In Sec. l VII we calculate dynamical properties from the thermo- 
dynamic pressure and we close with some concluding remarks 



in Sec.lVTII 



II. THERMODYNAMIC FORMALISM 

Our starting point is the dynamical partition function, 
which weighs points in phase space by the local stretching 
factors A(r, v, t) for trajectory bundles starting at phase space 
points (r, v) and extending over a time t. The stretching and 
contraction factors characterize the behavior of an infinites- 
imal volume in phase space under the dynamics. Typically, 
this volume will grow in some directions and shrink in others. 
Then the local stretching factor is the factor by which the pro- 
jection of the volume onto its unstable (expanding) directions 
will increase over time t. Similarly, the contraction factor is 
the factor by which the projection of the same volume onto its 
stable (contracting) direction is decreased over time t. Then, 
the partition function is defined by 



Z{f3,t) = / dfi{r,v)[A{r,v,t)]'-'', 



(1) 



where a temperature-like parameter f3 is introduced to empha- 
size the similarity to ordinary statistical physics. The integra- 
tion is over an appropriate stationary measure. 

From Eq.Q the topological pressure, P{(3), is obtained as 



P{(3) = lim - In Z{f3,t). 

t— too t 



(2) 



On introducing the Laplace transform Z(/3, z) = 
C{Z{P,t)}, the calculation of P{/3) greatly simplifies. 
From the definition of the Laplace transform we see that 
Z{(3, z) only converges if z stays smaller than a partic- 
ular value, the radius of convergence, which is given by 
z = limt^oo 7 In ^(/3, t) = P{P)- Thus, the topological 
pressure is given by the leading singularity of Z{P, z), ||7|. 

In analogy to the standard procedures of statistical physics 
we can define a dynamical entropy function, h{(3), as the Leg- 
endre transform of the topological pressure, i.e. 
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HP) = p{(i) - p 



dP{P) 
8/3 ' 



(3) 
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For long times the local stretching factors are approxi- 
mately given by the exponent of the sum of positive Lya- 
punov exponents, A^, multiplied by time, A ~ exp(i J^t ^i)- 
Therefore, it can be shown that the entropy function defined 
by Eq.Q can be identified for special values of /3 with dy- 
namical properties. For P — 0, h{(3) equals the topological 
entropy, /itop, whereas for (3 — \ \l equals the KS entropy, 
hKS which equals the sum of positive Lyapunov exponents 
Y^t A^, see for instance iillll- 

For systems where trajectories can escape, escape rates can 
be extracted from the topological pressure. E.g., for /3 = 1 
the topological pressure P{j3) equals —7, where 7 is the es- 
cape rate while the relationship — 9P(/3)/9/3|/3=i = \i 
remains valid. The intersection point of P{P) with the /3-axis 
can be related to the partial Hausdorff dimension, i.e. the frac- 
tal dimension of a line across the stable manifold of the at- 
tractor (see |10] for details), while the intersection point of 
the tangent at P(l) with the /3-axis is associated with the par- 
tial information dimension UOlA IJ . 



III. FIELD DRIVEN LORENTZ GAS 

We study the thermodynamic formalism for the dilute, field 
driven Lorentz gas without escape. This model consists of two 
species of particles. Heavy, immobile particles of radius a are 
placed at random positions, while point particles of mass m 
and charge q move in between them. The interaction between 
light and heavy particles is modeled by elastic collisions and 
the heavy particles are not allowed to overlap. An external 
electric field introduces a force which accelerates the moving 
particles in the direction of the field and an isokinetic ther- 
mostat prevents the system from heating up indefinitely. This 
thermostat keeps the kinetic energy and thus the speed v of 
each moving particle constant Then, the equations of 

motions are given by 



-E 



(4) 



where a — qE ■ v/mv'^ = v^e/v, with e = qE/mv. This 
choice assures that the kinetic energy is kept constant. The 
electric field is taken to point in the x-direction. 

Utilizing a Lorentz-Boltzmann equation (LBE) approach 
Van Beijeren et aJ. have calculated the Lyapunov spectrum 
of field-driven dilute Lorentz gases in two and three spacial 
dimensions 1 1, 2]. They did so by extending the usual LBE so 
as to include the so-called radius of curvature (ROC) tensor, 
p, (see Sec.|^). The resulting extended Lorentz-Boltzmann 
equation (ELBE) for the probability distribution of r, v, and 
p reduces to the usual LBE by integrating over p. From the 
ELBE the KS entropy is obtained as a steady state ensemble 
average. 



IV. RADIUS OF CURVATURE 

In order to describe the local stretching factor entering the 
partition function we have to measure the separation of two 




(a) without field 



(b) witli field 



FIG. 1: Illustration of the separation of two nearby trajectories, f{t) 
and r* (t), for (a) no external field and (b) with an external field, see 
text for details. 



nearby trajectories in the course of time. Sinai has given ge- 
ometric arguments relating this to the ROC |5, 13, 14, 15^. 
Since the light particle moves in d dimensions, and in tan- 
gent space the direction of the flow is neither expanding nor 
contracting, the unstable manifold (the subspace of expanding 
directions on the energy shell in tangent space) is of dimen- 
sion d — 1 (provided the driving field is not too strong). It 
may be represented by the set of all infinitesimal displace- 
ments in position space, orthogonal to the direction of the 
flow; the stretching factor for large t may be identified with the 
expansion factor of an infinitesimal volume in this subspace. 
This expansion factor in turn may be expressed in terms of the 
ROC. 

Fig^ shows a schematic plot of the radius of curvature in 
two spacial dimensions (in configuration space) in the cases 
with and without an external field. We measure the separation 
in a plane perpendicular to the reference trajectory. Let the 
spacial difference between two trajectories be dr{t) and the 
difference in velocity 5v{t) after some time t. These vectors 
satisfy the differential equation 



dSfjt) 
dt 



= Sv{t). 



(5) 



The solution of this defines the ROC tensor through the rela- 
tionship 



Sf{t) = -p{t,Sr{0)) ■ Sv{t). 

V 



(6) 



In the absence of external fields, with Sv being constant, the 
ROC tensor simply is of the form p{t) = vtl + p(0). Es- 
pecially in d — 2, where p reduces to a scalar, p{t) simply 
is the distance of the two nearby trajectories to their mutual 
intersection point, as illustrated in Fig[2 Its sign is positive if 
the intersection is located in the past and negative if it is in the 
future. 

One sees that already under free streaming the dependence 
of the ROC on its initial conditions becomes relatively less 
important as time increases. This is even much more so, if 
the free streaming is interrupted by collisions. In that case the 
ROC just before and just after a collision are related in d = 2 
by 



a cos ( 



(7) 
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which may easily be generalized to higher dimensions 
Here p_ is the pre-coUisional ROC, p+ the ROC directly after 
the colhsion, and the angle between the velocity vector at 
the collision and the outward normal to the scatterer at the 
point of incidence. Now, the point to notice is that for dilute 
systems typically pZ^ is of the order of the inverse mean free 
path, which is much smaller than the second term on the right 
hand side of Q. Therefore, to leading order in the density 
of scatterers this term may be ignored. This implies that the 
initial value of p already gets washed out after one collision. 
If one cannot use this low density approximation, still a few 
collisions suffice to make the ROC independent of its initial 
value. In the sequel of this paper we will use low density 
approximation, so we set 



_^ = 2/acos( 



(8) 



Combining Eqs. (|5} and (|6j one finds that that the stretching 
factor over a time r now may be expressed in terms of the 
ROC tensor as 



A(r) — exp 



V I dtdet(p"i(t)) 





(9) 



From this the KS entropy may be obtained as 



h-KS 



lim — 

t — 'OO t 



dTdct(p-i(T)) -«(det(p-i))^ 



(10) 

The bracket (■ • • )stat denotes an average over a stationary 
non-equilibrium distribution on phase space. Its equivalence 
to a time average requires ergodicity of the motion of the mov- 
ing particle on the chaotic attractor About the validity of this 
even less is known than about ergodicity in equilibrium, in the 
absence of a driving field. In our calculations we will actu- 
ally make plausible assumptions about the time average rather 
than using the phase space average. 

For calculating the increase of the stretching factor mea- 
sured from directly after a collision to directly after the sub- 
sequent one, one has to integrate det(p~^(i)) over this time 
interval and insert the result into Eq. (|5Jl. For N uncorre- 
cted collisions taking place over time t, the corresponding 
stretching factor is given by the product of the N individual 
stretching factors. The KS entropy is calculated by dividing 
the logarithm of the stretching factor by t and then taking the 
long time limit |l|l[l3. 



V. THERMODYNAMIC FORMALISM FOR FIELD 
DRIVEN LORENTZ GAS 

Van Beijeren and Dorfman present calculations in d 
dimensions for the dilute random Lorentz gas without an ex- 
ternal force. There, the dynamical partition function is calcu- 
lated by assuming that subsequent collisions are completely 
uncorrected. Under this assumption all free times between 
subsequent collisions may be assumed to be distributed ac- 
cording to the same exponential function and all scattering 



angles also follow the same simple distribution. In the present 
case, where the moving particle still has constant speed, we 
will make the same assumptions, but now the free motion in 
between collisions is not along straight lines any more. 

For simplicity we from now on will restrict ourselves to the 
case of d = 2, although the generalization to higher dimen- 
sions is fairly straightforward. As in (l^ we divide up the 
time interval [0,t] into subintervals [tr-i,tj], with i,. the in- 
stant of the r-th collision of the moving particle within this 
time interval. Let the total number of collisions be M, then 
to = and tAA+i = t- The total stretching factor can be 
factorized into a product as 



A(t) 



A/" 
r=0 



(11) 



where A^ = A(rr) is given by Eq.(|9j, with = — tr- 
For obtaining an explicit expression for this we have to solve 
the differential equation describing the time evolution of the 
ROC, which is of the form Ql 



p — V + pe cos a + p — sm 

V 



(12) 



where 9 is the angle between the velocity vector v{t) and the 
external field. As initial condition we will use the low density 
approximation pr^ — a/2 cos^r, with 0^ the collision angle 
of the rth collision. From the equations of motion it follows 
that 9 obeys the differential equation 9 = —e sin 9 which has 
the solution 



9{t,9^ 



arccot 



sinh(ei + K) , k = In [ cot(6'+/2)] . 

(13) 

Here, 0+ — 9{0, 0+) is the angle between the external field 
and the velocity direction directly after the rth collision, com- 
pare Fig[nb). With this solution for 9{t) we rewrite Ea.( ll2t 
as 



p = V + peta,iih{et + k) + p"^ 2 ■ (14) 

V cosh (et + k) 

We find for the ROC 

I ^\ 11 cosh K sinh K — ep^ 

, , , 1/1 V / V cosh K+ep-i- sinh k 

pit) — — cosh et + K,) r 

e 1_l_o^n'^^/'^^■/-_L^-^''^ cosn K smn K— ep+ 



1 + sinhfet + r — , ^-r— 

^ ' V cosn K.+ep^ smn /t 

(15) 

Inserting this into Eq.(|9j we find, to second order in e. 



A?) 



exp 



A, 
A(i) 



— cos 6*; 



A(i)A(2) 

VTr 



Pt+ 

24 



'1-9 sin^ 91 



(16) 
(17) 

.(18) 



Note that the above equations may also be used for r = but 
with Pg and 9'^ given by the initial values of p respectively 9 



at t = 0. 



For obtaining the topological pressure one has to substitute 
the above results into Eg. (II 1> . raise this to the power (1 — /3) 



4 



and average over the configurations of scatterers. For simpli- 
fying this average it turns out to be useful to rearrange Ea. Jl 1> 
as 



A(t) = «r^A(2) 



where 



-1.0+) 



(-■+1)+ 



-A?). 



(19) 



(20) 



are defined by 

[m}{P,z,e+') - J d0+M{l3,z,0+',e+)f{9+i25) 

{mff]{p,z,-) = f d9+Mfi/3,z,e+)f{e+) (26) 



The reduction to the independent variables implied here fol- 
lows from Eq.(|8jl combined with the relationship 



2(j)r = 9+ -9- ±11 



(21) 



and 9^ may be expressed in terms of 9+_^ and r^-i through 
Eq.O. 

A further simplification can be made by passing to the 
Laplace transform z) of the dynamical partition func- 
tion. Since, in the limit of large t, the topological pressure 
equals the logarithm of the stretching factor per unit time, it 
may be identified as the rightmost singularity of this Laplace 
transform. It has to be real, as the stretching factor is real 
and positive definite. Assuming independence of all free flight 
times Tj and scattering angles one finds straightforwardly 
that z) may be obtained as 



where integration limits of the -integrations depend on the 
sign of TT in the delta function in Eq.(l23ll and are given in the 
Appendix. Further, 1 denotes the unit operator and [• • • 
the operator inverse. 

From the above equations the rightmost singularity of Z 
readily follows as the value of z for which the largest eigen- 
value of the operator 271 equals unity. For zero field this eigen- 
value problem is trivial: the leading eigenfunction is the unit 
function, the eigenvalue is obtained easily and the resulting 
topological pressure coincides with that found in |16]. For 
small non-zero field one has to proceed by expanding OT, the 
leading eigenfunction / and the leading eigenvalue ji in pow- 
ers of the field strength, i.e. 



f 



(27) 
(28) 
(29) 



Then $H/ — jif is solved by standard perturbation methods 
using a Fourier series expansion for /, i.e. 



oo oo oo 

Z{f}, z) = / d^^^ • • • / d^o+ MfiP. 0+^) = + E cos{n9+) + ^ 6« sin{n9+), (30) 

\r_n n— 1 n=l 



M=0 



M=0 



1-971 



P„tj(/3-^--) (22) 



with i = 0, 1, 2. 

The details of the solution of this eigenvalue equation are 
given in the Appendix^ Here we just give the resulting eigen- 
value p to second order in the field strength, 

^i{(3,z) - A**°)(/3,z) + eV2)(/3,z) 

r(i)r(^)r(2-/3) 



Here the operators / and M are defined as the Laplace 
transforms of the configurational averages of the appropriate 

powers of vTj^KyJ and T{Tr, 9r+i, 9r), respectively. Specifi- 
cally, one has 



2^r(f + i) (i. + z)2-/3 

, (/3-l)(/3-2)(/3-12) 

48(1^ + 2)2 



(31) 



M{P,^,0+',9+) 



tt/2 

-7r/2 



Now the Laplace transform of the dynamical partition func- 
tion is of the form 



Z{f3,z)^C{p.,z) [l~^l{|3,z)Y 



(32) 



X S 



1 



Q+' 



- 9{t, 9+)±Tr) 



X [T{r,9+\9+)]''^ , 

OC 



(23) 

1-/3 



Mfif3,z,9^) 

with 9{t, 9+) defined in Eq.O- The operators and OT/ 



where the additional prefactor C, originating from 
z, •, 9^^), contains no singularities. 
If we also want to obtain information about the contract- 
ing direction, e.g. the negative Lyapunov exponents, we need 
to know the local contraction factors. For this we may con- 
sider the time reversed motion. Note, that we will consider 
the time reversed motion on the attractor for the "forward in 
time motion" and not the attractor of the time reversed mo- 
tion. The contraction factors can be calculated by considering 
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contracting trajectory 
bundle 



expanding trajectory 
bundle 



FIG. 2: Illustration of expanding and contracting trajectory bundles 
for the two dimensional field driven random Lorentz gas. 



contracting trajectory bundles instead of expanding ones fl4ll . 
see Fig|2| We keep the same notation as for the forward in 
time motion, therefore colHsion r precedes colUsion r — 1 in 
the time reversed motion. Hence the boundary condition for 
p{t) in Ea.( ll4t is p(t) = /9_ = acos(/)/2, where p_ speci- 
fies the ROC directly before a collision (in the forward in time 
motion). The ROC still evolves in time according to Eq. (ll4> 
for < t < T. Like the expansion factors, one calculates the 
local contraction factors by using Eq.(|9}. To second order in e 
this yields 



A,7 
All) 



Pr- 



-A? 



Al^ 



exp 



■ cos 



24 



P (l + 3sin2 0+) 



(33) 
(34) 

(35) 



with pr_ = a COS 4>r/'2- In the limit of vanishing external 
field the contraction factor is just the inverse of the stretching 
factor. This is no longer true if we apply an external field. We 
see from Eqs. ( I18t and ( I35t that differences occur in the field 
dependent exponential. 

To obtain a topological pressure for the contracting direc- 
tions we have to solve a similar eigenvalue problem as for the 
expanding directions. Since for the field free case the inverse 
of the contraction factor equals the stretching factor, we solve 
the eigenvalue problem for the inverse of the contraction fac- 
tor. The resulting eigenvalue, which has to be set equal to 
unity again, to second order in field strength becomes 



i-^r(i)r(/^)r(2-/?) 

2/5r(| + i) (l. + z)2-/3 
2(/3~l)(/3-2)(/3-^ 



48(i/ 



(36) 



Details are given again in the AppendiXn 

In the following section we will discuss the resulting topo- 
logical pressure and related properties. 



VI. DYNAMICAL PROPERTIES 

As stated in the preceding section the topological pressure 
follows as the value of z for which z) = 1. To second 
order in e this leads to 



P(/3)=Po(/3) 



48(Po(/3) + i') 



(37) 



where 



2 Va 



r(2-/3)r(i)rte 



2/3r ( § + 1 



(38) 

is the field free value of the topological pressure. The dynam- 
ical entropy then follows from Eq.Q as 



h{f3) = P(P) - (3 



dp 



12-/32 



48(Po(/3) + 1^) 
/3(/3-l)(/3-12) dPoiP) 



8(3 



(39) 



where again ho{(3) is the field free value of the dynamical 
entropy. 

We can perform similar calculations for the contracting di- 
rection. Then we obtain the equivalent of the topological pres- 
sure and the dynamical entropy. Note that in the limit of van- 
ishing external field the topological pressure for the contract- 
ing direction equals the one for the expanding direction. We 
obtain 



p-(/3) = PoW+e 



h-{l3) = hoiP) 



48(Po(/3)+z^ 



(40) 



48(Po(/3)+z.) 



8-^2 



/?(/3-l)(/3-8) aPo(/3) 



pm 



df3 



(41) 



Fig|3lshows the topological pressure and the dynamical en- 
tropy as functions of the parameter f3. As expected for a sys- 
tem without escape, P(/3) vanishes for (3 — 1. We further see 
that P(/3) is a convex function. 

A number of interesting dynamical quantities are related to 
the dynamical entropy, h{(3), for special values of /?. The KS 
entropy, Hks^ is given by the dynamical entropy for (3=1. 
From Eq.(l39ll we obtain 



(42) 



aiy\\ 11 

For the negative Lyapunov exponent we yield, from Eq. J4U 

7 e2 



-;i-(l) = A- = -z.(l-7-ln(^)) 



(43) 
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As one should expect, these resuhs coincide with those 
of previous calculations QJ 01 based on the same assump- 
tions (basically, no correlations between subsequent colli- 
sions). But here the results are extended to general values 
of /3. 

New results follow for the topological entropy, which is 
given by the dynamical entropy for f3 = 0. From Eq.(l39> 
we obtain 



''top 





1 

2 


e2 

— ly 

4 


" a ' 


1 

2 


. a . 









(44) 



That is, we obtain the zero field limit results given in fl6ll with 
a correction which is quadratic in the field strength. 

The equivalend of the topological entropy for the contract- 
ing direction is obtained from Ea. (l41> for (3 = 0, i.e. 



h-{0) 



6 



(45) 



We further calculate the dynamical entropy as a function of 
E{l3) = -dP{f3)/dp, 111). E{P) equals the average of the 
logarithm of the local stretching factors, (InA)^, where the 
subscript /3 refers to the fact that in phase space initial points 
are weighted according to the stretching factors raised to the 
power (1— /3). For (3=1 this yields the average of the positive 
Lyapunov exponent. 

The maximum entropy, which is the topological entropy, is 
always found for /? = as can be seen from Eq.Q- The KS 
entropy is given by the value for j3 = 1 where h[E{(3)] = 
E{(3), again according to Eq.(|3ji- In Fig@]a plot of h[E[(3)] 
vs. E{f3) is given. But notice that the descending parts of 
h[E{(3]) correspond to values of /3 < 0. 

The dynamical entropy is related to a dimension spectrum, 
D{a), with a = E{P), through D{a) = h{a)/a O- FigE 
shows the dimension spectrum, D[E{(3)]. From this one can 
find the partial Hausdorff dimension, Dh, as the maximum of 
D[E{l3y\. For systems where trajectories cannot escape, Dh 
equals one. This is clear because the maximum of D[E{(3)] 
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FIG. 3: The topological pressure, P{f3), and the dynamical en- 
tropy, h{P), for small fields and small densities of scatterers, i.e. 
e = 0.001, n = 0.001, v = 1, a = 1. 




0.01 0.02 
E(P) 



FIG. 4: Dynamical entropy, h[E{P)], as a function of E{l3) = 
-dP{P)/dp for e = 0.001, n = 0.001, a = 1, and = 1. The 
lower panel is a close-up of the boxed region of the upper panel. 




0.1 0.2 0.3 0.4 0.5 
E(P) 

FIG. 5: Dimension spectrum for e = 0.001, n = 0.001, a = 1, and 
V = 1. The maximum occurs for /3 = 1. 



is obtained for P{f3) = 0, i.e. (3=1, for systems without 
escape. Note that possible values of partial dimensions are 
restricted to the interval [0, 1], see SecHIl Furthermore, the 
full Hausdorff dimension is D^'' < J2'i=i ^ti where the 
D^^ are the partial Hausdorff dimensions corresponding to 
all stable and unstable manifolds iTs^l . 

The above results for the dynamical entropy allow to 
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calculate, respectively approximate, another dimension, the 
Kaplan- Yorke dimension, D^y In general the Kaplan- Yorke 
dimension is given by Dky — n + ^i/\^n+i\ where 

n is given by the largest value for which J^^^i — 0- Thus 
for the two dimensional Lorentz gas with constant energy we 
have one exponent equal to 0, one positive, and one negative, 
i.e. we have Dky 2 + A+/| |. There are only three ex- 
ponents since the system is restricted to a three dimensional 
hyperplane by the thermostat. From the KS entropies we see 
that Dky « 3 - 3e'^/{8h"j^giy), where /i^cs "^^e KS en- 
tropy without external field. Since the exact full Hausdorff 
dimension and therefore also the full information dimension 
cannot be extracted from the partial dimensions we assume 
that the Kaplan- Yorke conjecture still holds. This is consistent 
with our present knowledge: because -D^'' < X]i=i ^ti = 
we still can have I?^" = Dky- Let us further mention that 
the dimensional loss due to the external field, as expressed by 
Dky, is rather small in the region where the above results 
hold, i.e. for small fields and low densities, see also [19}. 

The calculated dynamical properties allow for the extrac- 
tion of macroscopic transport coefficients. The diffusion co- 
efficient is given by P = — \in\^Q[v'^{X+{e) + X-{e))]e^'^ = 
(3/8)?;^/i/, 01^ which can also be expressed in terms of the 
Kaplan- Yorke dimension by 2? = — \im^^Q[v^h'j^g{DKY ~ 
3)]e-2,Cl. 

Some comment on the calculations of the topological pres- 
sure as a function of the temperature-like parameter f3 is in 
order here. The obtained results have to be taken with a 
pinch of salt. In 1 16] and |20] it is shown that for the ran- 
dom Lorentz gas the results obtained there are restricted to 
increasingly smaller neighborhoods of /3 « 1 for increasing 
system size. As can be seen from Eq.Q, for /5 = 1 all points 
in phase space are equally weighted. For [3 < 1, though, the 
partition function will be dominated by the largest stretching 
factors, which correspond to the most unstable trajectory bun- 
dles. That is, for /? < 1 stretching factors from regions in 
phase space with a high density of scatterers, and therefore 
large stretching factors, dominate. So, in the limit t — + cx) 
it is possible that Z{P, t) is dominated by trajectory bundles 
confined to a small part of phase space. In regions of high 
scatterer density however, subsequent scattering events can- 
not be regarded as independent any more and the distribution 
of free times between scatterers in these regions will be very 
different from that for the system as a whole. With increasing 
system size the effects become more pronounced because the 
probability of finding approximately trapping regions of high 
scatterer density increases with system size. 

VII. CONCLUSION 

In the present study we have calculated dynamical proper- 
ties for the field driven random Lorentz gas within the thermo- 



dynamic formalism. In the limits of vanishing external field 
or (3 approaching unity, our results are in perfect agreement 
with those of previous studies. 

From the topological pressure we extracted various quanti- 
ties, such as the KS entropy and the topological entropy. A 
dimension spectrum was obtained by calculating the dynami- 
cal entropy as a function of the variable E{0), defined as the 
derivative with respect to (3 of the topological pressure. 

Van Beijeren and Dorfman have calculated KS entropies 
and topological entropies for general dimensions d for the ran- 
dom Lorentz gas without an external field llT3l . Presently we 
work on the extension of the present study to higher dimen- 
sions. Subtleties exist because an external field complicates 
the analytic calculation of the determinant of the inverse of 
the ROC tensor for higher dimensions quite a bit. However, 
the stretching factors can also be calculated by looking at the 
time evolution of the deviations in velocity, see 11411 . 

We also extended our studies to systems with open bound- 
aries ll2lll . This allows for using the thermodynamic formal- 
ism to study the escape rate formalism jE2il . as well as dimen- 
sion spectra for systems with escape. In the limit of ^ 1 
comparison can be made again with previous results IUTII . 
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APPENDIX: CALCULATION OF THE DYNAMICAL 
PARTITION FUNCTION 

Starting with the dynamical partition function we can de- 
termine the topological pressure which is given as the leading 
singularity of the Laplace transform of the dynamical partition 
function. 

In order to calculate the Laplace transform of the dynami- 
cal partition function, Eq. (l22> . we need to calculate the eigen- 
functions and eigenvalues of Tl. According to Ea. (l27> we 
expand in powers of e, yielding 



8 



7r/2 



2 



--(0+'-0(T,0+)±7r 



dr 

-7r/2 

X <^ 1 - e-cos0+ + [l - 9sin6l+ + Scos^ 6*+] 



2z;r 



1-/3 



(A.l) 



The operator 9Jt has to be understood as acting on a function f{0^). In order to eliminate the (5-function we integrate over 0+ 
first by noticing that we can write S[g{9^)] — S{9^ — Oq) / g' [Oq) where 9^ is the root of g{9Q) and the prime denotes the 
derivative with respect to 9^ . Here we have 



g'iet) = i - y cos 0+ + COS(20+) 



with 



(A.2) 



9+ =9+' 



- 2(1} + er sin 9:^ - -^-^ sin(26'+) ± tt. 



(A.3) 



Then the operator IDt(/3, z, 6'+') becomes 



m{l3,z,9 



7r/2 

dr j d(f> v e~ 

-Tr/2 



2VT 



acos( 



1-/3 



2 1 + er COS ( 



2 ^ 1-/3 

1 -e^cos6l+ +6^^ [-2 + 6cos(26'+)] 



(A.4) 



This is expanded to second order and the operator 971 acting on f{9^) then reads 



mi9+) 



(9+') = [^("Vl {9+')+e\m^^'^f] {9+') + [e2OT(2)/l (0+') 



(A.5) 



with 



mi")f] (0+') 



oo Tr/2 



dr 



-7r/2 



-(i'+z)t 



2ur' 



1-/3 



cos^ 0/(00+ 



oo 71-/2 

OT^^V (6'+') = - J dr J dcpve 

-7r/2 
oo ir/2 



-(i^+z)t 



2ur 



1 1-/3 



cos'3 0-(l + /3)cos0o+/((?o+) 



dr j 10 V e 

-7r/2 



-(i'+2:)t 



2ur' 



1-/3 



COS'^ 



r2 1^3/32 + 13/3 + 8 1) 
24 ^ 8 



(A.6) 



(A.7) 



(A.8) 



For small e we can expand cos 9^ and cos(20^) which gives Since we expand only to second order in e it is sufficient to 



cos 61+ = cos[6i+' - 2(/)±7r + ersin(6i+' - 20±7r)] 

« -cos(6l+'-2(/))-ersin2(6'+'-2(/)) (A.9) 



and 



cos(26i+) « cos(26i+' - 40) 



(A. 10) 



expand cos 9^ to first order because it only enters in the e- 
term of Ea.( IA.5> . Accordingly, we only take the zeroth order 
term of cos{29q). Now the eigenvalues /i and eigenfunctions 
/ can be calculated by using standard perturbation theory and 
a Fourier series expansion for / as given in Eq.( I30> . To zeroth 



9 



order in e we find that /''''^(6'+) = constant, which we set 
equal to 1. For the eigenvalue we get 



(0) 



a 



2/3(z/ + z)2-/3r(f + 1 

Inserting these results into Eq. jA.7t we get 

^(1) = 

and /W^) = -^r^cose 



(A. 11) 



(A. 12) 
(A. 13) 



The procedure for the e^-terms is analogous. However, 
when calculating the eigenvalue /i^^^ one finds an additional 

contribution from OT'-*'-'/*'^'' (^'^')' which also has to be 
taken into account. This is because of the e-dependent term 
in Eq.( IA.9l l. Then to second order in e the eigenvalue is given 
by 



(/3-l)(/3-2)(/3-12) 



(A. 14) 



In principle the eigenfunction /f^) (6/+) can be calculated and 
will be proportional to cos(2^^+). However, for our results we 



do not need Z^^) (5/+) since the only terms entering in Ea. J3H 

are^(")/(o)(0+)and eV^'V^"^^^)- 

An analogous calculation for the contracting direction 
yields for the eigenvalues 



-,(1) 



7(2) _ 



r(2-/3)r(i)r(^) 

A f(A.15) 



2'3(!y + z)2-/5r(f + 1) 



-,io) 



(/3-l)(/3~2)(/3-8) 



(A. 16) 
(A.17) 



and for the eigenfunctions 



r"\o+) = 1 



/ 



) 



4(z/ + z) 



cosf 



(A.18) 
(A.19) 



where the bar is indicating the contracting direction. Again, 
for Ea. (l36> we only need the eigenfunctions up to first order 
in e. 
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